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We present an effective mass theory for SiGe/Si/SiGe quantum wells, with an emphasis on calcu- 
lating the valley splitting. The theory introduces a valley coupling parameter, v v , which encapsulates 
the physics of the quantum well interface. The new effective mass parameter is computed by means 
of a tight binding theory. The resulting formalism provides rather simple analytical results for sev- 
eral geometries of interest, including a finite square well, a quantum well in an electric field, and 
a modulation doped two-dimensional electron gas. Of particular importance is the problem of a 
quantum well in a magnetic field, grown on a miscut substrate. The latter may pose a numerical 
challenge for atomistic techniques like tight-binding, because of its two-dimensional nature. In the 
effective mass theory, however, the results are straightforward and analytical. We compare our 
effective mass results with those of the tight binding theory, obtaining excellent agreement. 
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I. INTRODUCTION 

Silicon heterostructures form the basis of numerous 
semiconductor technologies. To understand future sil- 
icon heterostructure devices that involve quantum ef- 
fects, one must consider the quantum states associated 
with the degenerate valleys in the conduction band struc- 
ture. In a strained quantum well, the valley degeneracy 
is two-fold. This degeneracy is lifted by the singular 
nature of the quantum well interface, with characteris- 
tic energy splittings of order 0.1-1 meV for the case of 
SiGe/Si/SiGe quantum wellsi 1 ' 2 ' 3 ' 4 ' 5 ' 6 Because this split- 
ting is comparable in size to the Zeeman splitting, the val- 
ley states can compete with spin states for prominence in 
quantum devices.- For emerging technologies like silicon 
spintronic s 8 ' 9 ' 10 ^ 1 ' 12 and quantum computing ) 13 ' 14 ' 15 ' 16 
it is therefore crucial to obtain a solid physical under- 
standing of the valley physics, and to develop a predictive 
theory for the valley states. 

The technological significance of valley states was rec- 
ognized long before the current interest in quantum 
devices! 17 ' 18 Early studies focused on bulk silicon, partic- 
ularly on the electronic states of shallow donors. In this 
case, valley splitting was known to originate from the sin- 
gular core of the donor potential, known as the "central 
cell."— In one theoretical approach, the many-body in- 
teractions associated with the central cell were projected 
onto a single-electron, effective mass (EM) framework, al- 
lowing the energy spectrum of the low-lying hydrogenic 
states to be computed from first principles ! 20 ' 21 Reason- 
able agreement with experiments was attained. How- 
ever, because of the complicated projection procedure, a 
general purpose EM theory was not obtained until more 
recently^ 

Related techniques were applied to the problem of val- 
ley splitting in a semiconductor heterostructure. There 
were several attempts to develop an EM theory, which 
is well suited for treating inhomogeneous conditions, in- 
cluding conduction band offsets and non-uniform elec- 



tric fields associated with modulation doping and top- 
gates. However, the resulting theories proved controver- 
sial, and many important valley splitting problems re- 
main unsolved. 

Atomistic approaches like tight-binding (TB) theory 
have recently emerged as important tools for calculating 
the valley splitting. 2 - These techniques have been applied 
to a range of heterostructure geometries ; 23 ' 24 ' 25 provid- 
ing crucial insights and predictions for experiments that 
that may be difficult to implement. For example, the val- 
ley splitting has been predicted to oscillate as a function 
of the quantum well widths These oscillations can be 
reduced, or even eliminated, by applying an electric field. 

Here, we develop an EM formalism, which corrobo- 
rates the atomistic results quite accurately, and which 
provides a simple physical explanation for the intriguing 
oscillations. Specifically, we show that the behavior oc- 
curs because of valley coupling interference between the 
top and bottom interfaces of the quantum well. In an 
electric field, the wavefunction is squeezed to one side of 
the quantum well, thereby eliminating the interference ef- 
fect. We also use the EM theory to move beyond simple 
one-dimensional (ID) geometries. For example, quantum 
wells grown on a miscut substrate represent an inherently 
2D problem. Here, we show that interference effects also 
play a crucial role for such miscut geometries, causing a 
strong suppression of the valley splitting at low magnetic 
fields. 

The paper is organized as follows. In Sec. UH we review 
the two predominant approaches to valley coupling. In 
Sec. IIII1 we describe an extension to the conventional ef- 
fective mass theory that provides a perturbative scheme 
to incorporate valley coupling. In Sec. IIV1 we use a tight 
binding theory to calculate the new input parameter for 
the EM theory - the valley coupling v v - as a function of 
the conduction band offset AE C . In Sec. El we apply the 
EM theory to a finite square well geometry. In Sec. IVI1 
we obtain an analytical solution for a quantum well in 
an external electric field. In Sec. I VIII we consider the 
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experimentally important problem of a two-dimensional 
electron gas (2DEG). In Sec. IVIIIi we study the valley 
splitting in a magnetic field, when the quantum well is 
misaligned with respect to the crystallographic axes. Fi- 
nally, in Sec. lIXl we summarize our results and conclude. 



II. EFFECTIVE MASS APPROACH 

Earlier EM approaches for the valley splitting in a het- 
erostructure include the "electric breakthrough" theory 
of Ohkawa and Uemura ) 26 i 27 ' 28 and the surface scatter- 
ing theory of Sham and Nakayamad A review of these 
theories and other related work is given in Ref. [18| . The 
Ohkawa-Uemura formalism leads to a multi-valley EM 
theory based on a two-band model involving the low- 
est two conduction bands at the T point of the Brillouin 
zone, \iT!(S)} and \T 15 (Z)). 3 ^ The resulting bulk disper- 
sion relation has a local maximum at the T point, and 
it exhibits two degenerate valleys, at roughly the correct 
positions in fc-space. The sharp confinement potential at 
the heterostructure interface produces a natural coupling 
between the two z valleys. For an infinite square well 
of width L, the Ohkawa-Uemura theory obtains a valley 
splitting of the form E v ~ sm(2koL) / L 3 , where k$i is the 
location of the valley minimum in the Brillouin zoned 
This result was later confirmed by TB theory. 23 i 24 The 
theory therefore captures the main qualitative aspects of 
the valley physics, with no additional input parameters 
besides those describing the bulk dispersion relation. In 
this sense, it is a first principles theory of valley splitting. 

However, some aspects of the Ohkawa-Uemura the- 
ory have been called into question. First, it has been 
criticized for its inaccurate description of the disper- 
sion relation near the bottom of the valleysd leading 
to quantitative errors. More importantly, the method 
relies on a closed EM description, which cannot easily 
incorporate microscopic details of the quantum well bar- 
rier. This contradicts the fact that the valley coupling 
arises from physics occuring within several angstroms of 
the interfaced Such distances are much smaller than 
any EM length scale, and cannot be accurately described 
within any EM theory. 

A physically appealing description of the heterostruc- 
ture interface has been put forward by Sham and 
Nakayama. These authors develop a theory in which the 
reflection, transmission, and valley-scattering of waves at 
a Si/Si02 interface is built directly into the Bloch func- 
tion basis states. Since the confinement is incorporated 
into the basis set, it does not also appear as an external 
potential in the envelope equation. Any additional po- 
tentials entering envelope equation (e.g., electrostatic po- 
tentials) are therefore smooth, and easily accommodated 
in an EM approach. The analytical results for the valley 
splitting are similar to Ohkawa-Uemura. Indeed, the two 
approaches have been shown to be closely relatedd 

The Sham-Nakayama theory has also been criticized. 
First, the theory is not self-contained - a single input pa- 



rameter a is introduced to characterize the microscopic 
width of the interface. Although Sham and Nakayama 
provide an estimate for a, the parameter is phcnomeno- 
logical. More importantly, the resulting EM theory is 
somewhat cumbersome, and cannot provide simple ana- 
lytical solutions for the heterostructure geometries con- 
sidered here. 

In this paper, we develop an EM theory which retains 
the desirable qualities of both the Ohkawa-Uemura and 
Sham-Nakayama approaches. We introduce a valley cou- 
pling parameter v v , which efficiently describes the valley 
coupling for any type of interface, and which enables sim- 
ple analytical results for the valley splitting that are in 
agreement with atomistic theories. 



III. EFFECTIVE MASS THEORY 

The EM theory of Kohn and Luttinger 17 provides an 
excellent description of electrons in a semiconductor ma- 
trix under the influence of a slowly varying confinement 
potential U(r)d Here, "slowly" is defined with respect 
to the crystalline unit cell of length a: 

U(r)/|W(r)| >a. (1) 

In practice, the EM approach is extremely robust, of- 
ten proving accurate well outside its range of validity. 
Indeed, the standard textbook descriptions of shallow 
donors and quantum wells are both based on an EM 
theory, 32 despite the singular nature of their confinement 
potentials. 

When the validity criterion (fT|) is not satisfied, it is a 
good idea to compare the EM results with microscopic or 
atomistic approaches, such as the TB theory of Sec. IIVI 
For GaAs quantum wells, the EM theory provides quan- 
titatively accurate results in most situations. The ap- 
proach only breaks down for very narrow quantum wells, 
or for high subband indicesd On the other hand, for 
indirect gap semiconductors like silicon, the EM theory 
must be extended if valley splitting becomes an impor- 
tant issue. A singular confinement potential causes valley 
coupling, and calls for a more sophisticated treatment. 
Here, we provide a discussion of both the general multi- 
valley EM approach, and of the valley coupling, which 
arises from a sharp quantum barrier. 

In the standard EM theory, the wavefunction for a 
conduction electron in bulk Si can be written as a 
sum of contributions from the six degenerate valleys. 
However, the lattice mismatch between Si and Ge in 
a Sii-^Ge^/Si/Sii-^Gej; quantum well causes tensile 
strain in the Si layer (assuming strain- relaxed SiGe). As 
a result, four of the six valleys rise in energy, while the 
two z valleys fall in energyd The strain splitting is on 
the order of 200 meV for the composition corresponding 
to x = 0.3 d Consequently, only z valleys play a role in 
typical low-temperature experiments. 

The EM wavefunction for strained silicon can then be 
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FIG. 1: (Color online.) Comparison of effective mass and 
tight binding results for the two lowest eigenstates in a 
Sio.7Geo.3/Si/Sio.7Geo.3 quantum well of width 9.5 nm. (Only 
half the eigenfunctions are shown.) The solutions correspond 
to the same orbital subband, but different valley states. Bot- 
tom: ground state. Top: excited state. Solid lines: effec- 
tive mass theory. Circles: tight-binding theory. Dashed line: 
quantum well boundary. 



expressed as 



j=±z 



(2) 



The constants ctj describe the relative phase between the 
two z valleys, with \a+ z \ = \a- z \ = 1/V2. The functions 
e » 2 Uk .(r) are Bloch functions, where k± 2 = ±kQZ are 
the conduction valley minima. We shall see that the en- 
velope functions Fj(r) are the same for the two z valleys. 

A central feature of the EM formalism, which is a con- 
sequence of Eq. JTJ), is that wavefunction separates neatly 
into atomic scale oscillations (the Bloch functions) and 
long wavelength modulations (the envelope function). 
This is apparent in Fig. [TJ where we show the EM and 
TB results for the two valley states corresponding to the 
lowest subband of a quantum well, treating the Bloch 
functions as described in Sec. HVl The two wavefunc- 
tions have the same envelope, but their fast oscillations 
are phase shifted by tt/2. Note that although the EM 
and TB approaches are fundamentally different (discrete 
vs. continuous), their wavefunction solutions are almost 
identical. Also note that the fast oscillations arising from 
Ukj (r) (not pictured in Fig. [I]) are commensurate with the 
crystal lattice, while the oscillations from e ikjZ are not, 
since fc is not at the Brillouin zone boundary. 

We can draw two main conclusions from the EM/TB 
comparison. First, the EM treatment correctly captures 
the essence of the subband physics, including both the 
long-wavelength and atomic scale features. Higher sub- 
bands (not pictured in Fig. [T]) are also described accu- 
rately. Second, the main difference between pairs of val- 
ley states does not occur in the envelope function, but in 



the fast oscillations. To leading order, the valley states 
are fully characterized by their valley composition vec- 
tors a — (a_z,aq- z ). Since the envelope function is in- 
dependent of the valley physics at this order, the a vec- 
tors may be obtained from first-order, degenerate per- 
turbation theory, by treating the valley coupling as a 
perturbation^ We now describe the perturbation the- 
ory, following the approach of Ref. [22|], where shallow 
donors in silicon were considered. We specifically avoid 
the question of what causes valley splitting, since this 
lies outside the scope of the EM theory. Instead, we in- 
troduce valley coupling phcnomcnologically, through the 
parameter v v described below. 

In the EM formulation of Fritzsche and Twose ; 37 ' 38 the 
strained silicon wavefunction, Eq. @, is determined from 
the equation 



= 



j=±z 



[flo + K(«)-e]Fj(«). 



Here, 



and 



H = T(z) + V QW (z)+V <t> (z), 



h 2 d ( 1 d 



2 dz \mi dz 



(3) 



(4) 



(5) 



is the one-dimensional kinetic energy operator. The lon- 
gitudinal effective mass m; is materials dependent, and 
varies from layer to layer in the heterostructure. How- 
ever, for Si-rich SiGe layers, m; depends only weakly on 
the composition. We therefore take mi ~ 0.92 mo to be 
a constant. Note that the transverse effective mass rrit 
does not appear in Eq. ([5]), since we initially consider only 
one-dimensional problems. In Sec. IVIII1 a more compli- 
cated, two-dimensional problem is studied, in which rrit 
appears. 

Three different potentials energies appear in Eqs. ([3]) 
and ([!]). Vqw(z) describes the conduction band offsets 
in the heterostructure. For a quantum well, Vqw(z) cor- 
responds to a pair of step functions. The electrostatic 
potential energy V^(z) describes any additional, slowly 
varying potential. Typically, V<j,(z) ~ —eEz is an elec- 
trostatic potential caused by modulation doping in the 
heterostructure or by external gates. 

The EM approximation breaks down near a singu- 
lar confining potential like Vqw(z), leading to a val- 
ley coupling. Very near the singularity, criterion |T]) is 
not satisfied, and it becomes impossible to fully separate 
the short-wavelength physics of the crystal matrix from 
the long-wavelength confinement of the excess electron. 
However, we may neatly capture the valley interaction in 
terms of an effective coupling potential V v (z), which van- 
ishes everywhere except within about an atomic length 
scale of the interface. The detailed form of such a cou- 
pling may be quite complicated^ However, because it is 
so strongly peaked, we may treat V v (z) as a S- function 
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over EM length scales. Indeed, the 5-function formula- 
tion arises naturally from some first principles theories of 
valley coupling at heterostructure interfaces.— We then 
have 



V v (z) = v v S(z - Zi), 



(6) 



where z, is the vertical position of the heterostructure 
interface. In Sec. IVIII1 we consider a case where Zj de- 
pends on the lateral position, x. However, in the other 
sections of the paper, we assume that Zj is constant. The 
valley interaction potential V v (z) plays a role analogous 
to the central cell potential for an electron near a shal- 
low donor ^ The valley coupling strength v v is a scalar 
quantity, which must be determined from experiments, 
or from atomistic methods like TB. 

At lowest order in the perturbation theory (zeroth or- 
der), we set V v {z) = in Eq. (|3|). Because of the fast 
oscillations associated with the exponential factors, the 
contributions from the two valleys are approximately de- 
coupled at this level, reducing Eq. (J3]) the conventional 
Kohn-Luttinger envelope equation: 



h 2 d 2 
2m; az £ 



F^(z) 

= e (0) F (0)( z ). 



(7) 



Here, the superscript ^ denotes an unperturbed eigen- 
state. Note that the effective mass is the same for both z 
valleys. The corresponding envelopes are therefore equiv- 
alent, and we shall drop the valley index. 

We now solve for a. and e in Eq. Q using first or- 
der perturbation theory. By replacing Fj(z) in Eq. ([3]) 
with its zeroth order approximation, left-multiplying by 
(z)e lk ' z , and integrating over z, we can express 
Eq. |3]) in matrix form: 



A* 



Ai 

e<°> + A 



OL- z 



Ot-z 

a+z 



(8) 



We have dropped small terms involving atomic scale os- 
cillations in the integrand. The perturbation terms are 
defined as follows: 



A 



Ai = 



V v (z) FV>(z) 



dz, 



z~ 2lk « z V v (z) F^(z) 



dz. 



(9) 
(10) 



Diagonalizing Eq. (|SJ) gives the first order energy eigen- 
values 



e ± = e(°)+A ±|Ai|, 
and the valley splitting 

£Z„=2|Ai|. 



(11) 



(12) 



In the valley basis (a~ z ,a+ z ), the eigenvectors corre- 
sponding to e± are given by 



a± = -L(e^ 2 ,±e-^ 2 ) 



(13) 



Here, ± refers to the even or odd valley combinations, 
and we have defined the phases 

e i6 = A 1 /\A 1 \. (14) 

For the case of a single interface at z = z,-, we obtain 

A = v v \F^( Zi )\ 2 , A 1 = v v e~ 2ik ^\F^(z i )\ 2 . (15) 

We see that it is the magnitude of the envelope function 
at the interface that determines the strength of the valley 
coupling. The a± vectors for this geometry are obtained 
from Eq. [fig]), with 



= 2fc Zi. 



(16) 



The valley coupling integrals in Eqs. © and (jTUJ) pro- 
vide a simple but economical characterization of the val- 
ley splitting. The coupling parameter v v , which we com- 
pute below, provides a means to incorporate important 
microscopic details about the interface. The utility of the 
present theory is demonstrated by the ease with which 
we obtain results in the following sections. The accuracy 
of the theory is demonstrated in terms of the agreement 
between the EM and TB techniques. 

We close this section by noting that several authors 
have criticized the Fritzsche-Twose formulation of the 
multi-valley EM theory ; 18 i 40 particularly because of the 
way intervalley kinetic energy terms are calculated. 2 ^ 
However, we emphasize that the present treatment in- 
cludes all first order corrections to the valley splitting. 
Since the intervalley kinetic energy is of order v 2 , it has 
a weaker contribution. The previous criticisms therefore 
do not apply here. 



IV. TIGHT-BINDING THEORY 

We now discuss a tight binding method for modeling 
heterostructures in silicon. Our main goal is to compare 
the solutions from such an atomistic technique with those 
of the EM theory, and to compute the valley coupling 
parameter v v , whose value cannot be determined within 
the present EM theory. We focus on the two-band TB 
model of Boykin et a/. ) 23 i 24 because of its simplicity. 

In the two-band model, the TB Hamiltonian includes 
nearest neighbor and next-nearest neighbor tunnel cou- 
plings, v and u, respectively. The values v = 0.683 eV 
and u = 0.612 eV are chosen such that (i) the bulk dis- 
persion relation e(k) has two valleys, centered at |A;| = 
ko = 0.82 (2ir/a), and (it) the curvature of e(k) at the 
bottom of a valley gives the correct longitudinal effective 
mass mi = 0.91 mo. The unit cell in this theory consists 
of two atoms, with separation a/4 along the [001] axis, 
where a — 5.431 A is the length of the silicon cubic unit 
cell. These parameters correspond to bulk silicon. A 
more sophisticated theory should take into account com- 
positional variations and strain conditions. However, for 
most situations of interest, the modified parameters differ 
only slightly from the bulk. 



5 



> 

CD 




Band offset, A E c (eV) 

FIG. 2: The valley coupling parameter v v and the quantum 
well "set-back" parameter S, denned in Eq. (|27[1 . as a function 
of the conduction band offset AE C . The results are obtained 
for a symmetric, finite square well by comparing the EM and 
TB theories. 



In addition to the tunnel couplings v and u, we can also 
include onsite parameters Xi , to provide a locally varying 
confinement potential. Both the conduction band offset 
Vqw(z) and the electrostatic potential V^z) can be ex- 
pressed as on-site terms. To avoid boundary errors, we 
choose a TB lattice much larger than the confined elec- 
tronic wavefunction. Diagonalization of the resulting TB 
Hamiltonian gives energy eigenvalues, from which we can 
calculate the valley splitting. 

The two-band TB theory describes silicon valley 
physics with a minimal number (2) of input parame- 
ters, which are both fixed by fitting to measured band 
structure parameters. More sophisticated techniques can 
provide numerical improvements, but they generally do 
not capture any new physics. In Ref. J23j, a comparison 
of the valley splitting between the two-band theory and 
a detailed many-band theory shows excellent qualitative 
agreement. The more accurate treatment gives results 
that are smaller by an approximately constant factor of 
25%. 

Some typical TB eigenstates for a finite square well 
are shown in Fig. [TJ The data points correspond to the 
squared TB amplitudes plotted at the atomic sites. Com- 
parison with the full EM wavefunctions requires knowl- 
edge of the Bloch functions in Eq. ©. However, to make 
contact with the TB results, we only need to evaluate the 
Bloch functions at the atomic sites. According to Eqs. ([2]) 
and (|2"rj)) (see Sec.fV] below), the low-lying valley pair of 
EM wavefunctions can be expressed as 

* ± (z) = i [u_ k0 (r)e-* k ° z ±u +k0 (r)e+^ z ] F(z), 

(17) 

where the ground-state alternates between + and — as a 
function of L. We can denote the two atoms in the TB 
unit cell as A and B, with the corresponding Bloch func- 
tions Uk (A) and Uk (B). By translational symmetry, we 



must have Uk (A) = ±Uk {B) . (Typically, we observe the 
— sign in our TB analyses.) The Bloch functions satisfy 
time-reversal symmetry, so that u k() {r) = u-k (r). Thus, 
defining Uk (A) = |u/c (A)|e* v and assuming proper nor- 
malization, we see that the squared TB amplitudes must 
fall on the curves 



_(z)| 2 = 2cos 2 (k z + tp)F 2 (z) (18) 
„{z)\ 2 = 2sm 2 (k Q z + tp)F 2 (z), (19) 



where all information about the Bloch functions is re- 
duced to the unimportant phase variable (p. Eqs. (|18|) 
and (flT)]) are plotted in Fig. [H setting tp = 0. We see 
that these analytical forms provide an excellent repre- 
sentation of the TB results. However, we emphasize that 
Eqs. (fT8"|) and (fT§|) are accurate only at the atomic sites. 
To describe the wavefunction between the atomic sites 
would require additional knowledge of the silicon Bloch 
functions. 41 ^ 2 ' 43 ! 44 

We can use the TB theory to determine the EM valley 
coupling parameter v v by comparing corresponding re- 
sults in the two theories. This is accomplished in Sec. [V] 
for the finite square well geometry, with results shown in 
Fig. EI 



V. FINITE SQUARE WELL 

We consider a symmetric square well with barrier in- 
terfaces at Zi — ±L/2, corresponding to a quantum well 
of width L. We assume that the two interfaces are equiv- 
alent, so the same valley coupling v v can be used on both 
sides. Using EM theory, the resulting valley splitting is 



E v = Av v F 2 (L/2)\cos(k L)\. 



(20) 



An analytical solution of the EM equations for the en- 
velope function of a finite square well can be obtained 
by matching wavefunction solutions at the interfaces, 32 
giving 



F(L/2) = 



where 



1 k w L + sm(k w L) 
kb 2k w cos 2 (k w L/ 2) 



-1/2 



(21) 



(22) 



The wavevector k w can be obtained numerically from the 
transcendental equation 



k ln — 



2miAE c 

h 2 



cos 2 (k w L/2). 



(23) 



Some typical results for the valley splitting as a function 
of the well width are shown in Fig. [3] 

An approximate solution for Eqs. (|2"0")) - (l2"3"l) can be ob- 
tained in the limit of a very deep or a very wide quantum 
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FIG. 3: (Color online.) Valley splitting in a finite square 
well. We compare effective mass results (solid curve) with 
tight binding results (circles), as a function of the quantum 
well widths L and L' , respectively (see text). The TB data 
points occur at integer multiples of the TB unit cell. We as- 
sume a Sio.7Geo.3/Si/Sio.7Geo.3 quantum well, corresponding 
to AE C = 160 meV. 



well. When n 2 h 2 /2miL 2 < AE C , we find that k w ~ tt/L 
for the first subband, leading to 



2t2 



F 2 {L/2) 



n 2 h 



miAE c L 3 ' 
The valley splitting is then given by 
4v v n 2 h 2 



E v 



mi 



AE r L 3 



cos(fcoL)|. 



(24) 



(25) 



This agrees with the dependence on well width obtained 
in Refs. [23| and [28|, for an infinite square well. 

Diagonalization of the perturbation Hamiltonian in 
Eq. |(5J) gives the ground (g) and excited (x) valley state 
a vectors 



a g = (l,-sign[cos(fcoL)])/\/2 (ground), 
ol x = (1, sign[cos(/coL)])/v / 2 (excited), 



(26) 



where we have defined sign [a;] = x/|x|. These results are 
equivalent to Eq. (H2J), up to an overall phase factor. We 
see that the oscillations in Fig. [3] correspond to alter- 
nating even and odd ground states [a + = (l,l)/\/2 and 
ol = (1,— l)/v2, respectively], as a function of of L. 
This alternating behavior has been observed previously 
in TB analyses^ 

We can use these results to obtain an estimate for the 
valley coupling parameter v v , by comparing the EM and 
TB results for the finite square well. For the EM case, we 
first solve Eq. (|2"5|) numerically, to obtain k w . We then 
solve Eq. (|2ip for F(L/2), finally obtaining the valley 
splitting from Eq. (|20p. For a given value of AE C , we 
fit Eq. (|20p to the numerical TB results, as a function 
of L, using v v as a fitting parameter. However, there 
is an ambiguity in relating the quantum well width, L, 
in the continuous EM theory to the discretized width, 
V = Na/2, in the TB theory, where N is the number of 
silicon TB unit cells in the quantum well. For example, 
is the interface located on the last atomic site in the 
quantum well, the first atomic site in the barrier region, 
or somewhere in between? We see that the two widths, 
L and L' , may differ on the scale of a single atomic layer, 



or a/4 = 1.36 A. To allow for this, we introduce a second 
fitting parameter S, as defined by 



L = L' + 5 a/4. 



(27) 



With these two fitting parameters, v v and 5, we obtain 
nearly perfect correspondence between the EM and TB 
theories, as shown in Fig. [3l 

In this manner, we can map out v v and 5 as a function 
of AE C , giving the results shown in Fig. [2] We find that 
v v {AE c ) is linear over its entire range, with 



v„ = 7.2 x IQ-^AEr 



(28) 



Here, v v is given in units of eV-m when AE C is expressed 
in units of eV. As described in Sec. IIV1 a many-band TB 
analysis obtains results for v v that are smaller by a factor 
of about 25%. 

From Fig. [21 we see that 6 ~ 1.1 forms a reason- 
able approximation over the typical experimental range, 
AE C ~ 50-200 meV. This corresponds to about one 
atomic layer, or half an atomic layer on either side of 
the quantum well. We can interpret 5 as the set-back 
distance for an effective scattering barrier which causes 
valley coupling. A similar interpretation was given for 
the parameter a in Ref. [29|. We see that this set-back 
distance increases for a shallow quantum well. 

Finally, we consider the asymptotic limits of our EM 
theory. In the limit AE C — > oo, corresponding to an 
infinite square well, Eq. (|25p becomes exact. Since E v 
does not vanish, we conclude that v v oc AE C in this limit. 
This is precisely the behavior observed in Fig. O 

In the limit AE C — > 0, corresponding to a shal- 
low square well, Eq. (f2"5"]) is not valid. Instead, we 
obtain F 2 (L/2) ~ AE c Lmi/h 2 , leading to E v ~ 
2AE c Lv v mi cos(koL)/h 2 . In this limit, we expect the val- 
ley splitting to vanish, but we can make no other predic- 
tions about v v . The numerical results, however, suggest 
that the linear dependence of v v (AE c ) extends smoothly 
to zero. 



VI. QUANTUM WELL IN AN ELECTRIC 
FIELD 

We now consider a quantum well in the presence of 
an electric field oriented in the growth direction. The 
geometry is shown in the inset of Fig. 0J For physically 
realistic fields, caused by modulation doping or electri- 
cal top-gates, the resulting electrostatic potential satis- 
fies the EM criterion, Eq. (JTJ) . We therefore proceed as in 
Sec. El using the numerical values for v v (AE c ), obtained 
for a symmetric square well. The electrostatic potential 
enters our analysis through the envelope equation l[5]). 
Although there are no exact solutions for the problem of 
a tilted square well, the approximate treatment described 
here provides analytic results that accurately reproduce 
the results of TB theory. 

We assume an electrostatic potential energy given by 
V<p(z) = —eEz, and a quantum well of width L and height 
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For simplicity, we have chosen a wavefunction tail that 
terminates abruptly. This approximation is satisfactory 
for a variational calculation, since the tail, which is expo- 
nentially suppressed, contributes very little to the energy 
expectation value. 

In our trial function, the parameter k accounts for the 
finite barrier height by allowing the wavefunction to ex- 
tend into the barrier region. The upward shift of the 
wavefunction in the presence of an electric field is given 
by z t . The solution (|29|) becomes exact for an infinite 
square well at zero field, suggesting that the trial func- 
tion will be most effective in this limit. We therefore 
define the small parameters 



FIG. 4: (Color online.) Effective mass results for the val- 
ley splitting in an electric field, as a function of the quantum 
well width. Five E fields are shown, from bottom to top: 
E = 0, 1.04, 2.08, 3.12, and 4.16 MV/m. The well width L is 
scaled by the TB unit cell size, a/2 = 2.716 A. Results are 
shown for a Sio.gGeo.2j/Si/Sio.sGeo. 2 quantum well, analogous 
to Fig. 1(a) of Ref. [231 ]. Inset: iJ-field geometry, with con- 
finement potential V{z) = Vqw(z) + V<l>{ z ) an d variational 
wavefunction F(z) given in Eq. (|29[1 . Here, Zb = Zt — 7T /k. 



2^2 



R a 7T 



2miL 2 AE c 



and 



eEL 
4A£ C 



(30) 



The energy expectation value, obtained from envelope 
equation |J5J), can be expressed in terms of the dimension- 
less variational parameters 



AE C . The top barrier of the well lies at z = 0. We 
consider the following variational envelope function: 



F(z) = 



-^mn[k(z-zt)] 




{zt-%<z< z t ) 
(otherwise). 



(29) 



?i = kL and 82 = kz tl 



giving 



(31) 



AE r 



4TT 2 y(e 2 - 7r/2)/0! + 2n[ 



){ - 47r 2 y(6>2 - tt/2)/0i + 2tt6»|/3 

I — 



+ 2 -7T) 3 ]/3 (0! + 



h <7T) 

9 2 > tt). 



(32) 



In the first case in Eq. (|32|) . the wavefunction extends 
past both barriers. For large enough electric fields (y > 
x), the wavefunction only extends past one barrier (see 
inset of Fig. 4). Minimization of e with respect to 9\ and 
6> 2 gives 



9i ~ 7r + y/x -y - \/x + y 
0! ~ 7r( V /a:)V3 



(8y 2 x) 



1/6 



(a: > y), (33) 
(»<!/), (34) 



where we have made use of the fact that 62 <C n. For the 
case x > y, we have also used 7r — 6± <C 7r. 

To compute the valley splitting, we use the results 
F 2 (-L) ~ 2(tt - 6»i - 6» 2 ) 2 /L and F 2 (0) ~ 26>|/L when 



x > y, and F 2 (0) ~ W\Q\I-kL when a; < y to obtain 



4ti„ e 



(x- y) + 



, — ifc L 



(x + y)| 



(x > y) 

(a; < y)- 
(35) 

These solutions are plotted in Fig. 2] as a function of the 
quantum well width for several different E fields. In the 
figure, we have evaluated Eq. (|3"5|) only at integer mul- 
tiples of the TB unit cell, to facilitate comparison with 
Fig. 1(a) in Ref. [23|. The correspondence between the 
EM and TB theories is quantitatively and qualitatively 
accurate, particularly for large well widths. 

In Fig. [U we see that the crossover to high field behav- 
ior corresponds to E v becoming independent of L. The 
crossover occurs when y > x or 



E > 



2TT 2 h 2 



mie 



L 3 ' 



(36) 
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At low fields, the valley splitting exhibits interference os- 
cillations, similar to Fig. [3] At high fields, the envelope 
function no longer penetrates the barrier on the bottom 
side of the quantum well. The top and bottom barri- 
ers then no longer produce interfering contributions to 
the valley splitting, causing the oscillations in E V (L) to 
cease. Since we have chosen a trial wavefunction with no 
tail, the crossover to high field behavior in Fig. 5] occurs 
abruptly. This is in contrast with the TB results where 
small oscillations can still be observed right above the 
crossover. 

We now study the asymptotic behaviors of Eq. ([3"5]) . 
In the zero field limit, y — ► 0, Eq. ([35]) correctly reduces 
to Eq. (|25|) for a symmetric square well. In the high field 
limit, we can ignore the bottom barrier entirely. For an 
infinite barrier, the problem is often analyzed using the 
Fang-Howard trial wavefunction.— To make contact with 
this approach, we shall now perform a modified Fang- 
Howard analysis, and demonstrate a correspondence be- 
tween the two results. 

The conventional Fang-Howard trial function does not 
penetrate the barrier region. However, for a finite barrier 
height AE C , it is the amplitude of the wavefunction at the 
interface which determines the valley splitting. We must 
therefore modify the Fang-Howard trial function to allow 
the wavefunction to penetrate the top barrier, similar 
to Eq. (J2U). An appropriately modified trial function is 
given by 



F( \ _ \ -^/b 3 /2(z - zt) exp[b(z - z t )/2] for z < z t 
y ' \ for z>z t ' 

(37) 

This is just the usual Fang-Howard variational function, 
shifted upward by z t . In the Fang-Howard approach, the 
wavefunction tail decays exponentially into the lower por- 
tion of the quantum well. This treatment is more physical 
than the abrupt termination assumed in Eq. (|29p . How- 
ever, the tail contribution to the variational calculation 
is insignificant, to leading order. Following the conven- 
tional Fang-Howard approach, but now using b and zt as 
variational parameters, we obtain 



2v v eE 
AE n " 



(38) 



As in the case of a symmetric square well, we see that 
E v cx AE^ 1 . We also find that Eq. ([38]) is equivalent to 
the high-field limit of Eq. ([33]) . Thus, the non-oscillating, 
high-field portions of the curves in Fig. [4] correspond to 
the "Fang-Howard limit" described by Eqs. (J3Z} and ([38]) . 
This agreement demonstrates that the variational form 
of Eq. ([2"9"]l is robust over the entire field range. 



VII. VALLEY SPLITTING IN A 2DEG 

In this section and the next, we consider problems of 
particular experimental importance. In both cases, con- 
ventional techniques like TB theory are somewhat cum- 
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FIG. 5: (Color online.) Typical 2DEG structure for a 
SiGe/Si/SiGe quantum well. (a) Heterostructure, left to 
right: strain-relaxed SiGe substrate and barrier, strained sili- 
con quantum well, strain-relaxed SiGe barrier, n + SiGe dop- 
ing layer, and SiGe spacer layer, (b) Conduction band profile, 
showing the envelope function F(z) and the Fermi energy Ef- 



bersome. However, the EM formalism leads to straight- 
forward solutions. We first consider a two-dimensional 
electron gas (2DEG) in a SiGe/Si/SiGe quantum well, 
which extends the single electron senarios we have stud- 
ied so far, by including many-body interactions. 

We consider the modulation-doped heterostructure 
shown in Fig. Ufa). In this structure, we assume that 
the charge is found only in the 2DEG and the doping 
layers. A more detailed analysis could also include back- 
ground charge and charge trapped at an interface. Be- 
cause the modulation doping field is large in a typical 
heterostructure, we will ignore the bottom barrier in our 
calculations. The conduction band profile is sketched 
in Fig. [5Jb). We treat many-body interactions using 
the Hartree approximation, as common for a quantum 
well. However, other many-body interactions can also be 
included ) i 45 using similar techniques to calculate the 
valley splitting. 

The modified Fang-Howard approach gives a reason- 
able approximation for an electron in a 2DEG. 18 How- 
ever, we must include the electron-electron interactions 
self-consistently. Within the Hartree approximation, 32 
the 2DEG charge density p(z) is given by 



p(z) = -enF 2 (z), 



(39) 



where n is the density of the 2DEG and F(z) is defined 
in Eq. (13"?]) . The charge density is related to the electro- 
static potential through the Poisson equation, 



d 2 c 



<H 



dz 2 



P 

6 



(40) 



where e is the dielectric constant of silicon. The bound- 
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ary conditions on Eq. (|40p are given by 



h 1 1 dz = (z — > — oo) 
(z = 0). 



(41) 



The second boundary condition anchors the energy of the 
confinement potential, Vq-w(z), at top of the quantum 
welli22 The corresponding electrostatic potential is given 
by 



4>h{z) 



en 
2be 



[\b 2 {z - z t f - 4b(z - zt) + 6]e 
-[b 2 z 2 + Abz, 



b(z—zi) 



(42) 



The variational parameters are determined by mini- 
mizing the total energy per electron, given by 



1 



(43) 



Here, V^{z) corresponds to the Hartree potential 
—e4>H(z). The factor of 1/2 in the Hartree term pre- 
vents overcounting of the interactions^ 2 - Note that, in 
contrast with calculations for Si/Si02 inversion layers, 
we do not include any contributions from image poten- 
tials in Eq. (|43|) . since the dielectric constants for Si and 
SiGe are nearly equal, and any other interfaces that could 
produce images are far away from the 2DEG. 
We evaluate Eq. (|4*3"|) , obtaining 



27.2 



H 2 b 



en 



8mi Abe V 8 



33 



2bz t 



AE C b 3 zf 
6 ' 



(44) 



where we have made use of the dimensionless small pa- 
rameters bz t and e 2 nj eAE c b. The first of these describes 
the shift of the wavefunction towards +z, while the sec- 
ond describes the relative magnitude of the electrostatic 
energy with respect to the band offset AE C . Minimiza- 
tion of e with respect to b and zt gives 



8h 2 



33miAE c 



and 



33e 2 nm/ 
8H 2 ' 



(45) 



We can now calculate the valley splitting for a 2DEG. 
Under the previous approximations, we obtain a very 
simple expression for the wavefunction at the top quan- 
tum well interface, 



F 2 (0) 



2eAE c 



leading to the valley splitting 



EL 



v v e 2 n 

TaeZ; 



(46) 



(47) 



[Note that this result is obtained for a perfectly smooth 
interface. As we show in Sec. I Villi substrate roughness 
can reduce this estimate considerably. However, the scal- 
ing dependence in Eq. (I47[) remains valid.] 

We can use our estimate for the valley coupling param- 
eter, Eq. I|28p. to obtain a quantitative prediction for the 
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FIG. 6: (Color online.) Tilted quantum well geometry, with 
crystallographic axes (x,y,z) and rotated axes (x',y',z r ), as- 
suming y = y . The effective tilt angle is and the atomic 
step size is s. 



valley splitting in a 2DEG. Expressing the valley split- 
ting in units of meV and the 2DEG density n in units of 
10 12 cm -2 , we find that 



E„ 



l.Un. 



(48) 



Here, we have used the low-temperature dielectric con- 
stant for silicon, e = 11.4eo. It is interesting to note that 
the barrier height AE C does not directly enter the final 
result. 



We can compare Eq. (|48j) with the corresponding, 
non-self-consistent calculation, by treating the system of 
2DEG and doping layer as a parallel plate capacitor with 
electric field E — en/e. The valley splitting for an elec- 
tron in such a field is given in Eq. (|38[) , with the result 



E v ~ 2.29 n 



(non-self-consistent) . 



(49) 



The factor of 2 difference with Eq. ([48|) arises because 
the electric field in a real 2DEG is not uniform, due to 
the presence of charge. The non-self-consistent procedure 
therefore uses an electric field that is too large, overall, 
and it overestimates the valley splitting. 

In Eq. (|4"5)l . the prefactor 1.14 can be compared with 
similar estimates for silicon inversion layers. Ohkawa and 
Uemura obtain a prefactor of 0.15 ) 26 i 27 while Sham and 
Nakayama obtain 0.334^ A recent experiment in a top- 
gated Si02/Si/Si02 heterostructure has obtained a much 
larger value of the valley splitting4£ However, several re- 
marks are in order. First, we do not specifically consider 
top-gated structures here, so Eq. (|48|) does not directly 
apply to the latter experiment. Second, we point out that 
depletion and image charges, which were not considered 
here, play a more significant role in an inversion layer 
than a quantum well geometry. Finally, we note that a 
more accurate, many-band estimate for v v would reduce 
the prefactor in Eq. (|4"5| by about 25%, as discussed in 
Sec. HV] 
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VIII. TILTED QUANTUM WELL IN A 
MAGNETIC FIELD 



first express the envelope function equation in the tilted 
basis {x',y',z'), giving 



In high-mobility Si/SiGe devices, substrates are often 
intentionally tilted away from [001] by up to several de- 
grees, to promote uniform epitaxial growth under strain 
conditions, and to help reduce step bunching. Addi- 
tional roughening occurs during the growth of strained 
heterostructures. (Recent advances in nanomembrane 
technology may help overcome this difficulty. 48 ) Typical 
devices therefore contain atomic steps at their surfaces 
and interfaces, associated with global or local tilting. At 
the EM level, one might expect to ignore such small, 
atomic-scale steps, and to work in the locally tilted basis 
(x',y',z') shown in Fig. [S] However, when valley cou- 
pling is taken into account, the problem becomes two- 
dimensional, since the tilted surface is misaligned with 
respect to the crystallographic axes. In such high dimen- 
sional geometries, there is an obvious scaling difficulty 
for atomistic theories. However, the EM theory of valley 
splitting has a definite, practical advantage. The fol- 
lowing discussion expands upon our previous analysis in 
Ref. gjj. 

We consider a tilted quantum well in a magnetic field. 
In the low- field limit, the electronic wavefunction covers 
many steps, leading to interference effects in Eq. (fT0| . 
and a near-total suppression of the valley splitting. To 
see this, we note that the vertical position of the interface, 
Zi, becomes a function of the lateral position, x. The 
phase factor in Eq. (fT0|) is therefore not a constant, in 
contrast with the case of a flat interface. In Fig. we 
assume that a single step corresponds to a change of one 
atomic plane, or a/4 = 1.36 A. From Eq. (fT0|) . the phase 
shift between neighboring steps is given by 2fco(a/4) = 
0.82 7r. So the steps are nearly 180° out of phase, and the 
interference is severe. 

At low magnetic fields, the experimental data are con- 
sistent with this picture of suppressed valley splitting. 6 
(The experiments also suggest a non-vanishing zero-field 
extrapolation, which we shall investigate elsewhere.—) 
Some of the characteristic features of the magnetic field 
dependence can be explained in terms of the lateral con- 
finement of the electronic wavefunction. Magnetic con- 
finement over the length scale I b = \Jh/\eB\ reduces the 
number of steps that contribute to the valley splitting 
integral. The interference effects and the suppression of 
the valley splitting are similarly reduced. To get a sense 
of the scales involved, we note that a wavefunction of 
width 41b will cover 16(?s/a) tawd steps, where •& is the 
local tilt angle in Fig. [6j (Here, we have assumed a wave- 
function of diameter I^J (r 2 ) = 41b, obtained using the 
solutions described below.) For a typical ■& = 2° miscutj 6 
the wavefunction covers about 26 steps when B = 1 T, 
and 13 steps when B = 4T. The confinement provided by 
electrostatic top-gates can also enhance the valley split- 
ting by reducing the step coverage^ 

We can gain insight into the magnetic field dependence 
of the valley splitting by considering a simple model. We 



1 / d 

J2n i~ih— + eA n (r') 

t— 1 2m„ ^ fix 



.n=l 



d_ 

dxL 



(50) 



-Vqw(z') 



F{r') =eF(r'), 



where mi = TO2 = m t — 0.19 mo and 7713 = m; = 
0.92 mo, are the transverse and lateral effective masses, 
respectively. We note that the anisotropic effective 
masses are defined with respect to the crystallographic 
axes (x, y, z), not the growth axes. The full effective mass 
tensor in Eq. (|5"0|) should therefore include off-diagonal 
terms. In particular, we should have a term proportional 
to m~l ~ ^(m^ 1 — mj ), where we have taken d <C n/2. 
For a 2° miscut, however, we find that m t /m xz = 0.028, 
so the off-diagonal term is much smaller than the diag- 
onal terms. If desired, the off-diagonal corrections could 
be included, perturbatively. Here, we have considered 
only the leading order (diagonal) mass terms. 

In Eq. (fBTJ)) . the magnetic field is introduced through 
the vector potential. We consider the symmetric gauge, 
with A(r') = (—y',x',0)B/2. We also assume an ap- 
proximate form for the quantum well potential Vqw(z'), 
which is smoothly tilted (i.e., not step-like). This will be 
adequate for our simple estimate. Separation of variables 
then leads to solutions of form F(r') — F xy (x' , y')F z (z'), 
where F z (z') is the quantum well wavefunction, stud- 
ied elsewhere in this paper, and F xy (x',y') is the lateral 
wavefunction given by22 



F(f(r',6>) 



nl 2 B 2W+ 1 {n+ \l\)\ 



(51) 



xe 



L (W 



21% 



Here, we use radial coordinates, defined as (x',y') — 
(r' cos 6', r' cos 6'), while n — 0,1,2,... are the ra- 
dial quantum numbers (the Landau level indices), I = 
0, ±1, ±2, . . . are the azimuthal quantum numbers, and 
Ln l)[ \x) are associated Laguerre polynomials^ The en- 
ergy eigenvalues for F X y l \r\ 9') are given by 

2n + l + \l\ + l 
&ni — ru^c, 



(52) 



where lo c = e\B\/mt is the cyclotron frequency. Note 
that we have ignored spin physics here. 

To take an example, we now focus on the lowest 
Landau level, with n = I = 0. The behavior of the 
valley splitting in higher Landau levels is considered 
elsewhere M£l Eq. (JSTJ) now becomes 



F xy (x',y') = 



1 



(53) 
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Magnetic field, B (arb. units) 



FIG. 7: (Color online.) Valley splitting on a tilted quantum 
well, in a magnetic field. The lower (blue) curve shows the 
exponential suppression of the valley splitting, Eq. (|55j l , aris- 
ing from perfectly uniform steps. The upper (green) curve 
shows the linear behavior arising from the "plateau" model 
of disordered steps. Inset shows a localized step fluctuation 
(or "wiggle") of area A, as analyzed in the text. 



We assume a strong electric field, so that only one quan- 
tum well interface contributes to the valley splitting. For 
a smoothly tilted interface, the valley interaction poten- 
tial is given by 

V v (z') = v v 5(z'). (54) 



The perturbed valley splitting integral is therefore given 
by 

E v ~ \E v0 + Ae^Av v F 2 (0)F 2 y (r' ) sin(fc a/4)|. (56) 

Here, we have expressed the unperturbed result of 
Eq. (|55|) as E v0 , and we have assumed the amplitude of 
the envelope function is approximately constant across 
the wiggle at position r' = r' . We also note that the two 
terms in Eq. (|56[) enter the valley splitting integral with 
a phase difference <fr that depends on r' . Let us approx- 
imate F xy (r' ) ~ l/2wl% and sin(/coa/4) ~ 1. Then the 
fluctuation contribution is 

9 A 

«„i?(0)-3-, (57) 
nl B 

which has a linear dependence on the magnetic field. 
Therefore, at low fields, the fluctuation contribution 
dominates over the estimate obtained for a uniform tilt, 
Eq. (|55p. More realistic fluctuation models would in- 
clude a distribution of fluctuations, with contributions 
that partially cancel out due to interference effects. How- 
ever, the assumption of one dominant fluctuation loop 
(the "plateau" model^), leads to the linear dependence 
shown in Fig. which is consistent with experiments. 
(To obtain correspondence with Ref. [49[, we note that 
the "excess area" in that paper corresponds to 2A in our 
notation.) 



The valley splitting is then given by 



IX. CONCLUSIONS 



7T/| 



2 -r' 2 /Vl 



= 2v v F^0)e~ 2 ^ 1 ^ 2 



e 2ikoX '®6(z') dx'dy'dz' 
(55) 

-x' sin -d+z' cos -d ~ 



where we have used the fact that z 
— x'fi, along the z' = interface. 

The preceding results are obtained for perfectly uni- 
form steps at a quantum well interface, which we approx- 
imate by a smooth, uniform tilt. The resulting magnetic 
field dependence of the valley splitting, first reported in 
Ref. [49|, is shown in Fig. [71 The interference effect, aris- 
ing from the interfacial tilt, drives the valley splitting to 
zero at small fields, as consistent with experimental ob- 
servations. However, the exponential suppression of E v 
in Eq. (|55|) is an anomalous feature caused by the absence 
of disorder. If we consider more realistic step geometries, 
including disorder in the step widths and profiles^ the 
valley splitting will be enhanced by orders of magnitude, 
as confirmed by simulationsj 49 i 54 

We can obtain an estimate for the valley splitting en- 
hancement due to fluctuations by considering a single 
step wiggle, as shown in the inset of Fig. [7J The correc- 
tion to Eq. ((55]) is computed by noting that the area of 
the left step increases by A, while the right step decreases 
by A. We have shown that the phase difference between 
neighboring steps in the valley splitting integral is k^a/2. 



In this paper, we have developed an effective mass the- 
ory for the valley splitting of a strained SiGe/Si/SiGe 
quantum well. To compute the valley splitting of a per- 
fect quantum well with no steps, one needs two input 
parameters that describe the location and curvature of 
the band minimum as well a the valley coupling con- 
stant v v that captures the relevant microscopic details of 
the interfaces. These parameters must be obtained from 
a more microscopic theory, or from experiments. It is 
worth noting that, unlike bulk properties, the interface 
can vary from system to system, so that v v should be 
determined for each case. In this work, we have used 
a simple tight binding theory to compute func- 
tion of the conduction band offset AE C , assuming a sharp 
heterostructure interface. The results provide excellent 
agreement between the effective mass and tight binding 
theories. 

Once v v is known, we may apply the effective mass 
theory to a range of important problems. Here, we have 
considered the finite square well, with and without an 
applied electric field. We have also performed a self- 
consistent analysis of a 2DEG, using the Hartree ap- 
proximation. Excellent agreement between effective mass 
and atomistic theories confirms our main conclusion, that 
the valley splitting in this system can be fully explained 
through a single coupling constant, v v . 
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The effective mass theory is particularly useful for two 
or three-dimensional geometries, which cannot be easily 
treated in atomistic theories, due to scaling constraints. 
Here, we have applied the effective mass formalism to the 
inherently 2D problem of valley splitting in a magnetic 
field for a quantum well grown on a miscut substrate. We 
find that interference effects strongly suppress the valley 
splitting at low magnetic fields. However this suppression 
is reduced by introducing a small amount of disorder into 
the step-like geometry of the quantum well interface, as 
consistent with experimental evidence and simulations. 

N.B. In the final stages of preparation of this 
manuscript, we became aware of Ref. 55], which presents 



some similar results to those reported here. 
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